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Abstract. We present results for lattice QCD in the limit of infinite gauge coupling, 

^__( obtained from a worm-type Monte Carlo algorithm on a discrete spatial lattice but 

^^ with continuous Euclidean time. This is obtained by sending both the anisotropy 

^^ parameter 7^ ~ a/at and the number of time-slices N^ to infinity, keeping the ratio 

' r-l 1^ I^T — oT fixed. The obvious gain is that no continuum extrapolation N^ — > 00 has 

l~~5 to be carried out. Moreover, the algorithm is faster and the sign problem disappears. 

00 We compare our computations with those on discrete lattices. We determine the phase 

diagram as a function of temperature and baryon chemical potential. 

-)— > 
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D 1. Introduction 

The determination of the QCD phase diagram, in particular the location of the critical 

J> point, is an important, long standing problem, requiring non-perturbative methods. 

^ In lattice QCD, several approaches have been developed to investigate the phase 

10 transition from the hadronic matter to the quark gluon plasma, but all of them are 

limited to fiB/T<l, with hb the baryon chemical potential. The reason for this is 
f — . 
(^ the notorious sign problem, which arises because the fermion determinant for finite 

^~~' fi becomes complex, and importance sampling is no longer applicable. In QCD, the 

!| sign problem is severe, since the relative fluctuations of the complex phase factor grow 

exponentially with the lattice volume. In the strong coupling limit of lattice QCD 

(SC-QCD) discussed below, the sign problem is however mild and the full (/U,T) phase 

diagram can be obtained. 

2. Strong Coupling Lattice QCD 

In SC-QCD, the gauge coupling is sent to infinity and hence the coefficient of the 
plaquette term (3 = G/g"^ is sent to zero. Hence the Yang Mills part F^^F^i, is absent. 
Then, the gauge fields in the covariant derivative can be integrated out analytically. 
However, as a consequence of the strong coupling limit, the lattice spacing becomes 
very large, and no continuum limit can be performed. The degrees of freedom in SC- 
QCD live on a crystal. We consider the SC-limit for staggered fermions as studied in 
[1], which after gauge link integration and Grassman integration of staggered fermion 
fields yields a partition function that describes a system of confined, colorless, discrete 
degrees of freedom: 
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• Mesonic degrees of freedom k[i{x) G {0, ... 3} (non-oriented meson hoppings called 
dimers) and n{x) G {0, ...3} (mesonic sites called monomers), which obey the 
Grassmann constraint n{x) + J2u=±o ±d ^/i(^) = 3 at every lattice site x G A^^ x A^^; 

• Baryonic degrees of freedom, which form oriented baryon loops. These loops are 
self-avoiding and do not touch the mesonic degrees of freedom. 

The weights for temporal meson or baryon hoppings contain the anisotropy parameter 
7, needed in order to vary the temperature aT ~ 7^/iV^ continuously (see below). Here, 
we will restrict to the chiral limit, rUg = 0. In that case monomers are absent. 

3. Continuum Limit and Continuous Time Worm Algorithm 

Continuous time (CT) algorithms are now widely used in quantum Monte Carlo (see 
e.g. [2, 3]), but to our knowledge have not yet been applied to quantum field theories. 
The continuum limit in Euclidean time which we are interested in is the limit 

A^^ -> oo, 7 -)■ oo, I'^/N^ fixed (1) 

as 7^/A^r represents the temperature aT. Designing an algorithm that operates in 
this limit will have several advantages: There is no need to perform the continuum 
extrapolation A^^ — ?> oo, which allows to estimate critical temperatures more precisely, 
with a faster algorithm. Moreover, ambiguities arising from the functional dependence 
of observables on the anisotropy parameter (esp. non-monotonic behaviour as in Fig. 2 
left) will be circumvented. Also in the baryonic part of the partition function great 
simplifications occur: Baryons become static in the CT-limit, hence the sign problem is 
completely absent. The partition function can be written in terms of vertices at which 
spatial meson hoppings occur: 

Z{l,N^)^ E n vT^^\--^^^ n exp (-3a(x)/i,/T) (2) 

{k,cT} x(^Vm x(^Vb 

This relation becomes exact in the CT-limit since spatial dimers with multiplicity 2 or 
3 are suppressed by powers of 7 and are hence absent at 7 — )■ cxd, see Fig. 1 right: as 
the temporal lattice spacing at ~ a/7^ — ^ 0, multiple spatial dimers become resolved 
into single dimers. The overall number of spatial dimers remains finite in the CT-limit, 
as the sum over 0(7^) sites compensates the 1/7^ suppression. Temporal dimers can 
be arranged in chains of alternating 3-dimers and 0-dimers, which we denote by dashed 
lines, and 2-dimers and 1-dimers, which we denote by solid lines (see Fig. 1). The crucial 
observation is that the weight of these chains in the partition function is independent of 
their length, as the weight of each 3-dimer cancels that of the 0-dimer and likewise the 
weight of 2-dimers cancels that of 1-dimers. Hence, the weight of a configuration will 
only depend on the kind and number of vertices at which spatial hoppings are attached 
to solid/dashed lines, not on their position. For SC-QCD, there are two kinds of vertices, 
"L"-vertices of weight vl = 7~^, where dashed and solid lines join, and "T"-vertices of 
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weight vt = 2'j~^ /\/3, where a sohd hne emits a spatial dimer. The exponents ni{x) 
and nT{x) in Eq. (2) denote the number of T-vertices and L-vertices at spatial position 
X. In contrast to meson hoppings, spatial baryon hoppings are suppressed in the CT- 
limit by factors 7"^. Hence, baryons are static in continuous time. Positive (negative) 
oriented baryons are (dis)favored by a factor exp(±3/i/T) over meson lines. The sign 
problem has completely vanished! 

An important key step towards the CT algorithm is that spatial dimers are 
distributed uniformly in time as seen in Eq. (2). The lengths of dashed or solid intervals 
(which are related to the number of L- and T-vertices) are, according to a Poisson 
process, exponentially distributed: 



P(A/3) = exp(-AA/3), A/3 G [0, /3 = 1/aT] 



(3) 



with A the "decay constant" for spatial dimer emissions. Due to the presence of baryons 
A is space-time dependent: A = dM{x,t)/4:, where duix^t) is the number of mesonic 
neighbors at a given coordinate. Non-trivial meson correlations arise from the entropy 
of the various configurations. Likewise, baryonic interactions are due to the modification 
they induce on the meson bath, and thus also arise from entropy. 

The CT algorithm is a Worm-type algorithm, similar to the directed path algorithm 
introduced for SC-QCD in [?]. The updating rules are outlined in Fig. 1 and will be 
explained in detail in a forthcoming publication. 
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Figure 1. Left: updating rules for the continuous time algorithm. Right: illustrative 
2-dim. configurations (time flows to the right), in discrete time (top) and continuous 
time (bottom). Note how the latter lacks multiple spatial dimers and has only static 
baryons lines. 
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Figure 2. Left: continuous-time limit of the chiral susceptibility x at aT = 1.8 in the 
purely mesonic system U(3), exhibiting typical non-monotonic behavior in l/Nr- Note 
that the spatial suppression of double and/or triple spatial dimers has no effect on the 
CT-limit. Right: The SC-QCD phase diagram obtained on an Nt — 4 lattice [6] and 
in continuous time. 



4. Results on the SC-QCD Phase Diagram 

In SC-QCD at low temperatures, chiral symmetry, i. e. the Ua{^) symmetry of the one- 
flavor staggered action, is spontaneously broken according to Ul{1) x Ur{1) -^ Uv{l), 
and becomes restored at some critical temperature Tc. 

To obtain the SC-QCD phase diagram in the chiral limit mg = 0, we have measured 
the chiral phase transition temperature as a function of fi, and the nuclear transition 
with the baryon density. We were able to locate the tricritical point and also flnd a 
re-entrance already predicted by mean fleld analysis, due to the fact that the entropy 
decreases in the high-density phase as the lattice becomes saturated with baryons. Our 
new results eliminate systematic errors affecting previous flndings based on mean fleld 
approximations [5] or Monte Carlo for flxed A^,- [6]. 
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